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ABSTRACT 

We present results from Keck/NIRSPEC and Magellan/MMIRS follow-up spectroscopy of Lya 
emitters (LAEs) at z = 2.2 identified in our Subaru narrowband survey. We successfully detect Ha 
emission from seven LAEs, and perform a detailed analysis of six LAEs free from AGN activity, two 
out of which, CDFS-3865 and COSMOS-30679, have [O n] and [O ill] line detections. They are the 
first [O n] -detected LAEs at high-z, and their [Om]/[Oii] ratios and i?23-indices provide the first 
simultaneous determinations of ionization parameter and oxygen abundance for LAEs. CDFS-3865 
has a very high ionization parameter (qi on — 2.4T.g'g x 10 8 cms _1 ) and a low oxygen abundance 
(12 + log(0/H) = 7.84^q'25) in contrast with moderate values of other high-z galaxies such as Lyman- 
break galaxies (LBGs). COSMOS-30679 also possesses a relatively high ionization parameter (qi on = 
(3.4-15) x 10 7 cms _1 ) and a low oxygen abundance (12 + log(0/H) = 8.08-8.67), and falls below 
the mass-metallicity relation of z ~ 2 LBGs. Similarly, a low metallicity of 12 + log(0/H) < 8.4 is 
independently indicated for typical LAEs from a composite spectrum and the [Nil] /Ho index. Such 
high ionization parameters and low oxygen abundances can be found in local star-forming galaxies, but 
this extreme local population occupies only < 0.3% of the SDSS spectroscopic galaxy sample with a 
number density > 10 times smaller than that of LAEs. With their high ionization parameters and low 
oxygen abundances, LAEs would represent an early stage of galaxy formation dominated by massive 
stars in compact star-forming regions. Thus, low-mass galaxies such as LAEs produce ionizing photons 
efficiently probably with a high escape fraction achieved by density-bounded Hn regions, which would 
significantly contribute to cosmic reionization at z > 6. More than half of our LAEs have a larger 
EW(Lya) than expected from an instantaneous burst model at a given EW(Ha), suggesting that the 
majority of LAEs rather experience a more extended and continuous star-formation history. 

Subject headings: galaxies: evolution 



1. INTRODUCTION 



Lya emitters (LAEs), galaxies commonly observed at 
high redshifts with strong Lya emission, are considered 
to be low-mass, young galaxies as suggested from their 
small sizes, faint continua, and low masses inferred 
from spectral energy distribution (SED) fitting (e.g., 
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Vcnemans et al.ll2005t IGawiser et aT1l2006l : iPirzkal et all 
2007t lOverzier et al.l 120081) . LAEs are therefore likely 



to represent galaxies in the early st a ges o f galaxy 
evolution. More directly, iCowie et al.l (|2011l ) investi- 
gate rest-frame optical nebular lines of LAEs at low 
redshifts (z ~ 0.3). The authors find that a large 
portion of LAEs (75%) have equivalent width (EW) 
of Ha > 100 A, and that LAEs on average have lower 
metallicities and younger ages than the UV-continuum 
sample. These findings are consistent with the idea 
that LAEs are galaxies in early stages of galaxy forma- 
tion. At higher redshifts, LAEs are efficiently detected 
thanks to narrowband imaging techniques which have 
enrich ed our know l edge a bout the younger u n iverse 
(e.g.. ICowie fc Hul Il998t iMalhotra fc Rhoadsl 120021: 
IQuchi et al.l 120031: IGawiser et al. 2006; Shimasak u et al.l 
20061: IGawiser et all l2007t iGronwall et al.l 1200 



Ouchi et all l2008t iNilsson et all 120091: iGuaita et all 



2010; iFinkelstein et all 12011b iNakaiima et all 1201 



However, more direct evidence supporting the idea that 
high redshift LAEs are in an early evolutionary phase of 
formation is still needed. 

The gas-phase metallicity is a key property of galaxies, 
since it is a record of their star-formation histories. This 
physical quantity is relatively easily constrained with line 
ratios of nebular lines at rest-frame optical wavelengths 
(e.g., iNagao et al.l 120061 ). Another key quantity is the 
ionization parameter, defined as the ratio of the mean 
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ionizing photon flux to the mean hydrogen atom den- 
sity. Since the excitation of the Hn region is sensitive to 
the age distribution of the exciting stars, the ionization 
para meter provides a ro ugh estimate of age of a galaxy 
(e.g., IDopita et al.l I2006h . A large ionization parameter 
is observed from a galaxy dominated with massive stars, 
which is a sign that the galaxy is in an early stage of 
galaxy formation. Determining the ionization param- 
eter requires two emission lines of different ionization 
stages of a same elemen t, such as [On] and [Oin] (e.g., 
IKewlev fc Dopital[2002l ). 

However, well-defined samples of LAEs are generally 
located at very high redshifts (3 < z < 7), where 
rest-frame optical nebular line s are redshifted i n to in- 
frar ed wavelengths. Recen tly, iFinkclste in et al.l ([201 If ) 
and [Nakajima et al.l (|2012T) instead take advantage of 
LAEs at moderately high redshifts (z ~ 2), where the 
nebular lines up to Ha are observable from the ground 
in the near-infrared (NIR) windows. 

These studies extend the met allicity censuses to galax- 
ies with strong Lya emission. [Finkclste in et all ([201 If ) 
obtain NIR spectra for two z ~ 2.4 LAEs, and find 
that at least one LAE appears to be less chemically en- 
riched than z ~ 2 continuum-selecte d galaxies at sim- 
ilar stellar masses ([Erb et al.l l2006al ) . INaka iima et aLl 
(2012) estimate average line fluxes for [On]A3727 and 
Ha+[Nn]A6584 for z = 2.2 LAEs by stacking 1.18 and 
2.09 fim narrowband images for more than 100 LAEs, 
and place a firm lower-limit for the average metallicity 
of this population. Interestingly, the lower-limit is higher 
than expected for its stellar mass at the redshift. 

In contrast to the gas-phase metallicities, the ioniza- 
tion parameters of LAEs are unknown. Such measure- 
ments have only been obtained for bright galaxies such as 
Lyman-break g alaxies (LBGs) or strongly lensed galaxies 
at high-z (e.g . . iPettini et al.l 1200 H Mainline et all 120091 : 
iRichard et all 1201 If ). Although the current number of 
measurements is small, ionization parameters of such 
high-z galaxies are likely to be higher on average than 
those of local galaxies. A comparison of ionization pa- 
rameters for LBGs and LAEs will enable us to discuss 
whether or not LAEs are young galaxies at high red- 
shifts. 

In order to obtain reliable measurements of proper- 
ties including metallicity and ionization parameter of 
LAEs at high redshifts, we have carried out a large 
survey for z = 2.2 LAEs, using our custom narrow- 
band filter NB387 with Subaru/Suprime-Cam. At z = 
2.2, important nebular lines such as [On]A3727, H/3, 
[O III] AA5007, 4959, Ha, [Nn]A6584, are observable from 
the ground. Initial results were based on using three 
narrowba nds to detect Lya, [O III , and H a over the same 
vol umes (INakaiima et al.l l2012t see also iLee et all 120121 
and lLv et al.ll20lTl for the two NIR narrowbands and the 
NewHa Survey). 

In this paper, we present results from NIR spec- 
troscopy. We used Keck/NIRSPEC and Magel- 
lan/MMIRS spectrographs, and successfully detected 
Ha emission from seven LAEs. The number is dou- 
ble the previous number of high-z LA Es with NIR 
spectra (two from [ Finkclst ein et all 1201 ll and two from 
iMcLinden et alj2011l) . and allows us to begin to examine 
statistical variation of rest-frame optical spectroscopic 
properties of LAEs at z ~ 2. Our first NIR spectroscopic 



re sult discusses the kinem atics of LAEs and is presented 
in lHashimoto et all (|2012[ ). As a companion study, this 
paper presents the ionization and chemical properties of 
LAEs based on multiple nebular lines. Remarkably, we 
detected [O n] and [O in] lines from two LAEs , which 
provide ionization parameter estimates for LAEs for the 
first time. In addition, the oxygen lines allow us to deter- 
mine oxygen abundances that compliment the previous 
spectroscopic constraints on metallicity of LAEs from the 
[Nil] /Ha index. 

We also investigate Lya and Ha hydrogen lines 
for LAEs. A comparison of observed Lya/Ha ratios 
with the Case B r ecombination value (Lya/Ha = 8.7; 
iBrocklehurstl 1 1 9711 ) provides important insights into the 
physical mechanisms causing the strong Lya emission in 
LAEs. Furthermore, we study Lya and Ha equivalent 
widths which probe star-formation history, stellar age, 
and metallicity. 

This paper is organized as follows. We describe the 
data in ^5] The detection and measurement of emis- 
sion lines in the NIR spectroscopy is summarized in Sj3[ 
In Sj4[ we derive properties of LAEs including estimates 
of ionization parameter, metallicity, and SFR from the 
rest-frame optical nebular lines. We also check for the 
presence of active galactic nuclei (AGN) in the LAEs. 
In fJSJ we compare LAEs with other galaxies in terms 
of their ionization state, metallicity, and SFR. We then 
discuss the implications. We also discuss the physical 
properties inferred from Lya and Ha emission. We con- 
clude the paper in fj6] with a summary. Throug hout 
this p aper, magnitudes are given in the AB system (jOkd 
|1974f ) , and we assume a standard ACDM cosmology with 
(n m ,n A ,H ) = (O^OJ^Okms^Mpc- 1 ). 

2. NIR SPECTROSCOPIC DATA 
2.1. Sample Construction 

We carried out NB 387 imaging observa tion with 
Subaru/Suprime-Cam dMivazaki et al.l I2002D on 2009 
December 14-16 and 19-20 to search for z = 2.2 LAEs. A 
total of ~ 1.5 square degrees are covered, with pointings 
in the following five fields: the Subaru/XMM- N ewton 
Deep Survey (S XDS) field (iFurusawa et all [20081) . the 
COSMOS field dScov lie et all 120071) . theChandra Deep 
field South fCDFS: IGiacconi et al l [200l. the Hawa ii 
Hubble Deep Field North fHHDFN; iCapak et al.l l200l, 
and the SSA22 field (e.g., ISteidel et all 120001 ). The first 
results of th e NB387 survey i n the SXDS field have been 
presented in INaka iima et ahl (j2012f ). For the other fields, 
we select z = 2.2 LAE candida t es in the same manner 
as presented in INakaiima et al.l ([20121 ) , and a summary 
of the full ~ 1.5 square degrees survey will be presented 
elsewhere (K. Nakajima et al. in preparation). Briefly, 
objects which are bright in NB387 compared to the U 
and B bands are selected as LAE candidates. Interlop- 
ers such as [O n] emitters, Civ emitters, or active galac- 
tic nuclei (AGNs) are removed from the catalog by using 
UV, X-ray and/or radio data. 

Seven candidates were selected for NIRSPEC follow- 
up. We selected COSMOS-13636, COSMOS-30679, and 
COSMOS-43982, since their Lya spectra had been ob- 
tained with Magellan/MagE (M. Rauch et al. in prepa- 
ration). COSMOS-08501, HHDFN-18325, and HHDFN- 
18431 were selected because of their large Lya equiva- 
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TABLE 1 

Summary of the LAE sample with NIR spectroscopy 



Object 


R.A. 


Decl. 


NB387 C/-NB387 B-NB387 EW(Lyo) 


F(Lya) 

(5) 


EXPTIME< 6 ) Instr. 




(i) 


<i) 


(2) (3) (3) (4) 


(J) (H) (K) ( 7, 



COSMOS-08501 


10:01:16.80 


+02:05:36.26 


23.94 


1.41 + 0.09 


2.05 + 0.10 


255 ± 26 


24.0 + 1.1 




3600 


N 


COSMOS-13636 


09:59:59.38 


+02:08:38.36 


23.53 


0.91 + 0.05 


1.03 + 0.06 


73 + 5 


32.4 + 1.4 




5400 


N 


COSMOS-30679 


10:00:29.81 


+02:18:49.00 


23.63 


0.86 + 0.05 


0.55 ± 0.06 


34 + 3 


19.1 ± 1.4 


5400 


7200 6300 


N 


COSMOS-30679<t) 






23.73 


1.20 + 0.07 


1.16 + 0.07 


87 + 7 


23.2 + 1.8 








COSMOS-43982 


09:59:54.39 


+02:26:29.96 


23.83 


1.14 + 0.08 


1.30 + 0.07 


105 + 8 


30.3 + 1.3 




3600 


N 


HHDFN-18325 


12:36:23.36 


+62:06:05.10 


21.84 


1.26 + 0.02 


1.50 + 0.02 


122 + 2 


125. ± 1.2 




3600 


N 


HHDFN-18431 


12:36:25.62 


+62:05:37.43 


23.20 


1.15 + 0.03 


1.71 + 0.04 


156 + 6 


37.8 + 0.7 




3600 


N 


CDFS-3865 


03:32:32.31 


-28:00:52.20 


22.29 


1.42 ± 1.23 


1.01 + 0.42 


64 + 29 


84.0 + 13.9 


5100 


10800 


N/M 


CDFS-6482 


03:32:49.34 


-27:59:52.35 


23.26 


1.47 + 1.92 


0.96 + 0.66 


75 + 52 


41.9 + 22.0 




10800 


M 


SSA22-8043 


22:17:47.33 


+00:08:28.36 


24.51 


0.66 + 0.06 


0.52 + 0.06 


28 + 4 


6.3 + 1.1 




10800 


M 



Note. — (1) Coordinates are in J2000. (2) NB387 aperture magnitude. The diameter of the aperture is 2''0, and the typical error is 0.04 mag. (3) 
u*(U)— NB387 and B — NB3S7 colors and their la errors calculated from aperture magnitudes. (4) Rest-frame EW of Lya emission line in units of A 
calculated from B — NB387 color and rcdshift of Ha. For HHDFN objects whose Ha are not detected, we assume z — 2.18, which corresponds to the 
wavelength of the peak of the NB387 transmission curve. (5) Flux of Lya emission line in units of 10 -17 ergs -1 cm" 2 calculated from EW(Lya) and 
B band total magnitude. (6) Exposure time with NIRSPEC J, H, and K bands in units of second. For the CDFS and SSA22 objects, their H and K 
bands spectra were taken with MMIRS using the HK grism. (7) Instrument used for the spectroscopy. "N" stands for NIRSPEC, and "M" for MMIRS. 
For CDFS-3865, its J band spectrum was obtained with NIRSPEC, and H and K bands spectra with MMIRS. (f) Values obtained after removing the 
contribution from an adjacent object ( S|4.2I >. 



lent widths and fluxes measured from the NB387 imag- 
ing data. CDFS-3865, whose rest-frame optical spectrum 
has in part already been taken with MMIRS, was also 
observed with NIRSPEC for its [On] detection. Tabled] 
summarizes the details of the sample. 

For the MMIRS observation, we observed CDFS and 
SSA22 fields with one mask each. Details of the obser- 
vation and data reduct ion procedures are presented in 
iHashimoto et al.l {2012). Briefly, we found three LAEs 
in which we identified Ha emission. In this paper, we 
make use of the spectra of the three LAEs; CDFS-3865, 
CDFS-6482, and SSA22-8043 (Table [J). Figure ED shows 
the distributions of the LAEs presented in this paper on 
the U— NB387 versus NB387 color magnitude diagram. 

2.2. NIRSPEC Observation 

The observations were carried out on 9-10 February 
2011. Both nights were photometric. We observed in 
low-resolution mode, with a slit width of 0'.'76 and a 
slit length of 42'{3- We observed all six candidate with 
Nirspec-6 filter (hereafter referred to as K band) . In ad- 
dition, we observed COSMOS-30679 with Nirspec-5 (H 
band) and Nirspec-3 (J band), and CDFS-3865 with J 
band. The wavelength ranges of K , H, and J bands are 
1.88-2.31 /xm, 1.47-1.76 /im, and 1.15-1.36 /im, respec- 
tively for our configurations. The resolution in the K 
band is R ~ 1500. Exposure times for each object are 
given in Table [TJ 

Our science targets were acquired using the invisible 
object acquisition procedures. In this mode, an align- 
ment star brighter than ifvega ~ 18 was placed in the 
slit simultaneously with a science target. We first ac- 
quired the star at the center of the slit, then we nodded 
the telescope so that both the science target and the star 
were on the slit with the same distance from the slit cen- 
ter. For CDFS-3865 and HHDFN-18431, since we could 
not find any stars to place on the slit along with them, we 
first acquired the nearest star to the target (at distances 
of 68" for CDFS-3865 and 39" for HHDFN-18431), then 

11 The effective slit length is ~ 38", since we did not place 
objects in < 2" from the edges of the slit. 



nodded the telescope with an offset, which was calcu- 
lated from the WCS difference between the star and the 
science target. For HHDFN-18325, although a star was 
found at a distance of 32", which is smaller than the slit 
length, half of the observation was done without the star 
in order to avoid the persistence caused by the previous 
target alignment star. During an exposure, we manually 
guided the slit using the slit-viewing camera; if the star 
started to drift out of the slit, we manually moved the 
slit so that the star (and thus the invisible science target) 
stayed in the slit throughout each exposure. 

Standard stars, which were selected from the Hipparcos 
catalog to have AOV spectral type and similar airmass to 
the science targets, were observed at the beginnings and 
the ends of the nights. The calibration data were taken 
in an ABBA position pattern, while the science data in 
AB position pattern. 

2.3. NIRSPEC Data Reduction 

We used the Keck IRAF-based WMKONSPEC pack- 
ag43 to reduce the data. We also used IRAF scripts 
that were originally written for reducing near-infrared 
multi-object spectroscopic data fr om Subaru/MOIRCS 
(MCSMDP; lYoshikawa etHW^OWt) . 

The data reduction process included bias subtraction, 
flat fielding, distortion correction, wavelength calibra- 
tion, cosmic ray rejection, (A-B) sky subtraction, resid- 
ual background subtraction, image shifting, and stack- 
ing. We used the bias and flat data that were obtained in 
the same night as the science data. We performed wave- 
length calibration usi ng OH-lines by comp aring them 
with an OH-line list (jRousselot et al.l 120001) . For cos- 
mic r ay rejection, we used LA. COSMIC ( van Dokkuml 
120011 ). The sky and OH-lines of an A-position image were 
roughly removed by subtracting an average B-position 
image created from the previous and following images. 
We then removed the residual sky by subtracting a 9-th 
order polynomial fit in the spatial direction after mask- 
ing columns of positive and negative parts caused by 
the alignment star and the object. Individual images 

12 http : //www2 . keck . hawaii . edu/inst/nirspec/wmkonspec . html 
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Fig. 1.— Distribution of all objects in the u*(C/)-NB387 vs. NB387 plane detected in the COSMOS and SSA22 fields (left) and in the 
CDFS field (right). The red circles show the LAEs presented in this paper, and the black dots show the NB387-detected objects. For the 
purpose of display, objects whose u*(U)— NB387 colors exceed 2.4 are plotted at u* (£/)— NB387= 2.4. The horizontal solid line in each 
panel shows the selection threshold of u* (U)— NB387 and the blue curve indicates the 2<r photometric error in u* (U)— NB387 for objects 
with u*— NB387=0.05 (U— NB387=0.38), which is the average u*(U)— NB387 color of all the objects. The right j/-axis in each panel shows 
the rest-frame Lyc? equivalent width of z = 2.18 LAEs with the u*(V)— NB387 color corresponding to the left j/-axis. 



TABLE 2 

Fluxes of nebular lines determined from NIR Spectroscopy 



Object [On]A3727 H/3 [Om]A4959 [Oni]A5007 Ha [Nn]A6484 



COSMOS-08501 ■ • ■ ■ ■ • ■ • ■ ■ • ■ 1.91 ± 0.36 < 0.36 

COSMOS-13636 ■■■ ■■• ■■■ ■•■ 2.71 ± 0.38 < 0.38 

COSMOS-30679 1.04 ±0.26 ... (t) . . . (t) 4.11 ±0.41 3.11 ±0.27 < 0.27 

COSMOS-43982 ■■• ■■• ■•■ ■•■ 6.63 ± 0.58 4.14 ±0.58 

Composite^ ■■■ ■■■ ■■■ ■•■ 2.64 ±0.20 < 0.20 



CDFS-3865 2.53 ± 0.41 13.0 ±4.1 19.0 ± 3.5 53.6 ± 3.5 38.5 ± 2.4 < 2.4 

CDFS-6482 ■•• ••• ■■■ 23.8 ± 2.7 9.38 ± 1.71 < 1.71 

SSA22-8043 ■■■ ■■• ■■■ 11.0 ±4.2 14.7 ± 4.3 < 4.3 

Note. — Fluxes and their la errors are given in unit of 10~ 17 ergs~ 1 cm -2 . For lines with less than the 3a 
detection level, we list their la upper-limits. (*) Composite spectrum of the four objects f E|4. 11 . (f) These lines 
suffer badly from OH-lines subtraction errors. 



were shifted in spatial direction so that the object is in 
the same position both in A-position and B-position im- 
ages. The offset values were derived from the position 
differences of the alignment stars. When no star was 
observed simultaneously, we defined the offset values as 
the nod separation size. Finally, we stacked the position- 
matched individual images to create the two-dimensional 
(2D) spectra. 

We obtained flux solutions by comparing spectra of 
the A0V standard stars and a model, which was created 
by a stellar spectral sy nthesis prog ram (SPECT RUM; 
IGrav k CorbaiIv1[l99l based on the lKuruczl (fl99l 's at- 
mosphere models. The model spectra were then normal- 
ized so that their </, H , and K band magnitudes matched 
with those of standard stars, whose photometry was ob- 
tained from the Two Micron All Sky Survey (2MASS) 
All-Sky Point Source Catalog. The one-dimensional (ID) 
spectra were extracted by summing up 6-10 pixels in the 
spatial direction. The lengths of columns were deter- 
mined based on the seeing conditions; about two times 
the seeing size was used. We also confirmed that the 
sizes were large enough to detect most of the signal, and 
maximize the signal-to-noise (S/N) ratio. The ID and 
2D spectra are shown in Figure [5] We note in addition 
that the science targets and the standard stars used for 
the flux calibrations were observed and reduced with al- 
most the same conditions, thus effects of slit losses were 
also corrected in the procedures, since the standard spec- 



tra extracted in the procedure were normalized to have 
the total magnitude of that star. LAEs ty pically have 
half-l ight radii < f.5kpc (< 0'/18 at z = 2.2; IBond etaU 
I2009T ) . which is much smaller than the seeing size. There- 
fore, we cannot resolve LAEs with our observations, and 
assuming the PSF profile for the LAEs is reasonable. 

Emission line fluxes were measured by fitting a Gaus- 
sian profile to each line with the IRAF task splot. The 
sky noise level was estimated in the following manner; 
we prepared an aperture box, which had approximately 
twice the seeing size in spatial direction and twice the 
FWHM of the best fit Gaussian to the emission line in 
wavelength direction. We spread more than 100 aper- 
ture boxes around the emission line after masking pixels 
heavily contaminated by OH-lines, and measured their 
photon counts. We then fit the histogram of the counts 
with a Gaussian, and regarded its a as the ler fluctuation 
for the aperture used to measure the emission line. The 
line flues and their lcr errors are summarized in Tabled 

3. EMISSION LINE DETECTIONS 

3.1. Ha detection 

We detect significant Ha emission in K band 2D spec- 
tra for COSMOS-08501, COSMOS-13636, COSMOS- 
30679, COSMOS-43982, CDFS-3865, CDFS-6482, and 
SSA22-8043, but we do not identify any emission lines 
for the other two LAEs, HHDFN-18325 and HHDFN- 
18431. Both LAEs are expected to have strong Lya, 
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Fig. 2. — From left to right, J, if, and X band spectra of the LAEs. In each panel, the 2D and ID spectra are shown in the top 
and bottom, respectively. The vertical dotted lines in the ID spectra show the expected locations of nebular emission lines; [Oll]A3727, 
H/3, [O III] AA4959, 5007, [Nll]A6548, Ha, and [Nll]A6584. The gray shaded areas overlaid on each ID spectrum highlight spectral regions 
strongly affected by OH-airglow. For CDFS objects, their K band spectra at 2.115 £im and longer are not obtained due to out of sensitivity 
of MMIRS, and also shaded with gray. 
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Fig. 3.— BPT-diagram l IBaldwin et alj[l98lf) . The dotted and 
dashed curves show th e dividing line betwee n sta r-forming galaxies 
and AGNs defined by Kcwlcy ct al. (2001) and Kauffmann et al. 
(2003), respectively. CDFS-3865 is shown with the red circle, while 
the other LAEs are shown with the vertical solid lines, due to the 
lacks of [Om] and/or H/3. Since we do not identify [Nil] except for 
COSMOS-43982, upper-limits of [Nil] /Ha are given. The [Nil] /Ha 
ratio for COSMOS-43982 is given with a horizontal error bar. The 
red dot-dashed vertical line indicates the upper-limit of [Nil] /Ha 
for the composite spectrum of NIRSPEC-detecte d LAEs. The two 
triang les represent LAEs at z ~ 2.3 and 2.5 from Finkclstcin et al. 
I|2011l ). The gray small points are randomly selected objects from 
the SDSS spectroscopic sample. 

hence strong Ha. A possible reason why we do no de- 
tect lines for HHDFN-18431 is that no alignment star 
was observed simultaneously, and the object may have 
drifted significantly out of the slit. The same issue may 
explain the non-detection of Ha for HHDFN- 18325, since 
half of the observation was done without an alignment 
star f H2.2p . Since we have not yet obtained Lya spectra 
for the two objects, the non-detections may be alterna- 
tively due to the mis-selection of LAEs. Future optical 
spectroscopy is needed to resolve this issue. In the fol- 
lowing sections, we use the seven LAEs with reliable Ha 
detection. 



3.2. Oth 



er emission 



line detections 



For CDFS-3865, we additionally obtain [On]A3727, 
H/3, and [O iii]AA4959, 5007 emission lines. Combined 
with the Ha-detection, a full suite of prominent rest- 
frame optical nebular lines are thus obtained for this 
object. For COSMOS-30679, we obtain [On] and 
[Oin]A5007 emission lines, while we do not identify H/3 
and [Om]A4959, because they fall in a dense OH-line 
wavelength range. We thus cannot place meaningful 
upper limits on their fluxes. When fluxes for these 
lines are required for inferring physical properties, we 
assume a value of 0.28 for the [O iii]A4959/5007 ratio 
and calculate the H/3 flux based on the observed Ha 
flux, the intrins ic Ha /H/3 ratio a ssuming Case B recom- 
bination (2.86; lOsterbrockl 119891 ) . and the dust extinc- 
tion estimated from SED fitting. We emphasize that 
the two LAEs, CDFS-3865 and COSMOS-30679, are the 
first [O ii]-detected LAEs individually at high-z. The 
[Oii]-j-[Oili] lines provide the simultaneous determina- 
tions of ionization state and oxygen abundance (t j4.3|) . 

For CDFS-6482 and SSA22-8043, we detect 
[O in] A5007 while we do not identify H/3 and [O ill] A4959, 




Fig. 4. — K band composite spectrum of the NIRSPEC-detected 
objects. 

probably due to the contamination by OH-lines and 
the limited sensitivity of the instrument (More detailed 
descriptions about the MMIRS objects are provided by 
Hashimoto eTaT]|2"oT3 ). For COSMOS-43982, we clearly 
detect [Nn]A6584 emission line, which is not identified 
for the other LAEs. The [Nn]-detection suggests a 
non-negligible contribution of AGN to the emission 
lines for the object. Further investigations of AGN 
contamination are discussed in M4. 1 1 

4. PHYSICAL PROPERTIES OF LAES 
4.1. Removal of objects with AGN 

AGN, as well as star-formation, can produce large 
amounts of ionizing photons, and show strong Lya. Since 
our interest is star-forming LAEs, we need to remove 
LAEs with AGN activity from the sample. 

During the LAE sample selection, br ight LAEs with 
AGN signature have been removed (see INakaiim a et al.1 
[2011 . Obvious AGNs are therefore excluded from the 
sample. However, the procedure is not always perfect, 
and can miss objects with relatively weak AGN activ- 
ity. In order to assess possible contamination of AGN 
in our LAE sam ple, we plot our LAE s on a BPT di- 
agram (Figure [3] iBaldwin et al.l 119811 ) which is widely 
used to separate star forming galaxies from AGN. In 
Figure the underlying small gray points denote the 
spectro scopic objects fro m the Sloan Digital Sky Survey 
(SDSS: lYork et alJl2000f Pl The two curves shown in Fig- 
ure [3] are the empirical d emarcations (jKewlev et al.ll200ll 
Kauff mann et al.1 120031 ). Since only one LAE, CDFS- 
3865, has both [Om] and H/3, all the other LAEs are 
shown with vertical lines. Due to the non-detection of 
[N ii] , the vertical lines show upper-limits of [N n]/Ha ex- 
cept for COSMOS-43982, which has an individual [Nil] 
detection. Based on its relatively high [Nil] /Ha ra- 
tio, COSMOS-43982 may be contaminated by an AGN, 
though there remains a possibility that this LAE may 
be a star-formation dominated galaxy with a very small 
[O m] /H/3 ratio. For this analysis, we choose to be con- 
servative and regard the object as a candidate possess- 
ing AGN activity. The properties of COSMOS-43982 are 
considered and interpreted carefully in the following sec- 
tions. 

13 We use a part of the spectroscopic data taken 
from t he MPA-JHU DR7 release o f spectrum measurements: 
http : //www.mpa-garching.mpg. de/SDSS/DR7/ 
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TABLE 3 

Broadband Photometry of the LAEs 



COSMOS sample^ 


B 


V 


r' 


i' 


z' 


J 




K s 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


GUbMUb-UBoUl 


25.86 


25.91 


26.08 


25.88 


25.81 


98.45 




25.64 


99.99 


99.99 


99.99 


99.99 


COSMOS-13636 


24.43 


24.21 


24.35 


24.19 


24.24 


23.10 




23.43 


24.10 


23.75 


99.99 


99.99 


COSMOS-30679 


24.05 


23.12 


22.91 


22.46 


22.33 


21.15 




21.82 


22.12 


22.57 


99.99 


23.06 


COSMOS-30679<t) 


24.76 


23.82 


24.44 


24.09 


23.49 


22.31 




23.29 










COSMOS-43982 


25.00 


24.38 


24.48 


23.99 


23.73 


21.89 




21.62 


21.20 


21.02 


20.69 


20.75 


(limitmag) 


(29.13) 


(28.18) 


(28.33) 


(27.87) 


(26.89) 


(24.17) 




(24.84) 


(25.05) 


(24.25) 


(21.90) 


(20.63) 


(limitmag) (t) 


(28.76) 


(26.24) 


(25.79) 


(25.34) 


(24.85) 


(23.68) 




(24.63) 










CDFS sampled 


B 


V 


R 




z' 


J 


H 


K 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


CDLS-3865 


23.01 


22.94 


22.92 


23.14 


22.93 


22.73 


22.27 


22.38 


22.82 


22.82 


22.51 


23.00 


CDFS-6482 


23.93 


23.87 


23.78 


23.95 


23.67 


23.50 


23.36 


23.07 


22.88 


22.83 


22.34 


99.99 


(limitmag) 


(28.32) 


(27.85) 


(27.82) 


(26.14) 


(25.64) 


(24.57) 


(24.57) 


(23.97) 


(26.23) 


(25.68) 


(23.66) 


(23.43) 


SSA22 sample® 


B 


V 


R 


V 


z> 


J 




K 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


SSA22-8043 


24.80 


24.63 


25.54 


24.64 


24.59 


24.10 




21.79 


23.33 


23.14 


26.58 


21.15 


(limitmag) 


(27.85) 


(27.99) 


(28.03) 


(27.75) 


(27.16) 


(24.74) 




(22.78) 


(25.40) 


(23.64) 


(22.71) 


(20.95) 



NOTE. — Broadband photometry of the LAEs. All magnitudes are total magnitudes. 99.99 mag means no signal detected. Magnitudes in 
parentheses are lo - uncertainties adopted in SED fitting. 

U) BVr'i'z' data are obtained from Subaru/Suprimc-Cam, J data from UKIRT/WFCAM, K s data from CFHT/WIRCAM, and [3.6] - [8.0] data 
from Spitzcr/IRAC. All data are collected from the COSMOS Archive. 

< 2 ) BVRI data arc obtained from MPG/ESO 2.2m/WFI, z' data from CTIO/MOSAIC II, JK data from CTIO/IS PI, H data from NT T/SOFI, and 
[3.6] — [8.0] data from Spitzcr/IRAC. Optical and NIR data are collected from the MUSYC Public Data Release JGawiser et al.||2006l) . and Spitzer 
data from the SIMPLE Legacy II project. 

< 3 ) BVRi'z' data are obtained from Subaru/Suprimc-Cam, J data from K PNO/NEWFIRM K da ta from UKIRT/WFCAM, and [3.6] - [8.0] data 
from Spitzer /IRAC. Optical data are provided by T. Hayashino (see also iHayashino et al.l 120041) , J data are obtained from the NcwHq survey 
(Ly et al. 2 01l| ; [Lee"et alj|20l3) . K data are collected from the WFCAM data access page, and Spitzer data from the Spitzer Heritage Archive. 
^ Values obtained after removing the contribution from an adjacent object ( EI4.2B . The limiting magnitudes are calculated by adding the photometric 
errors and additional errors (residuals after subtracting GALFIT models) in quadrature. The limiting magnitudes are thus shallower than those 
given above (photometric only). 



For CDFS-3865, thanks to its relatively strong con- 
straint on [Nil] /Ha, we can see that the object is along 
the star forming sequence and has negligible AGN ac- 
tivity, similar to two other L AEs with spectroscopi c 
follow-up at similar redshifts (Finkclste in et al.l [2011). 
Although the other LAEs have no constraint on their 
y-axis values, AGN contributions to them are assumed 
to be minimal, since there are few AGNs in the range 
log([N n]/Ha) < —0.5, as suggested from the underlying 
SDSS galaxies. In order to check the assumption, we 
stack the K band spectra of the four NIRSPEC-detected 
LAEfLl with inverse- variance weights. The composite 
spectrum is shown in Figure SI and provides an average, 
deeper constraint on [Nil] /Ha ratio for the LAEs. The 
red dot-dashed line in Figure [3] shows the upper- limit of 
[Nil] /Ha ratio constrained by the composite spectrum. 
The line is further away from the area occupied with 
AGNs. Therefore, AGN activity is on average negligible 
for the LAEs. In the following sections, we regard all the 
LAEs except for COSMOS-43982 as galaxies dominated 
by star-formation. 

4.2. Stellar Mass and E(B-V) from SED fitting 

We perform SED fitting to broadband photometry to 
infer properties of the stellar populations. The ma- 
jority of images are collected from publicly available 
databases (see notes in Table [3]). The optical and NIR 
photometry is done with a 2'.'0 diame ter aperture using 
the d ouble- image mode of SExtractor ([Bertin fc Arnoutsl 
1996). The aperture magnitudes are then converted 

14 We do not use the MMIRS spectra due to their worse sensi- 
tivity and spectral resolution. 



into total magnitudes using aperture correction values, 
which are estimated from differences between aperture 
magnitudes and MAG_AUT0 values for point sources. For 
the IRAC i magin g, we use a 3'.'0 diameter aperture and 
lYan et all (|2005D 's aperture correction values. Table [3] 
summarizes the results. 

The procedure of the SED fitting is the same as that 
of lOno et al.l (|2010() . except for fixed redshifts, which 
are derived from Ha. Briefly, we u se the stellar popu- 
lation synthesis model GALAXEV (jBruzual &: Chariot! 
I2003T) for stellar SEDs, and include nebular emission 
(|Schaerer fc de Barrosl I2009D. A Salpeter initial mass 
function (IMF; iSalpeter] I1955D is assumed. We choose 
constant star-formation history and the stellar metallic- 
ity Z — 0.2 Zq . For dust extinctio n, we use Calzetti's 
extinction law (jCalzetti et al.1 12000]) on the a ssumption 
of E(B - V) gas = E(B - V)* as proposed by lErb et al.1 
(2006b), although this assumption is still debatable 
for high-z gal axies (e.g., iForster Schreiber et al.l [20091 : 
ILv et al] l2012t ). Intergalactic medium (IGM) attenua- 
tion is calculated using the prescription given bv lMadaul 
(|1995| ). These assumption s are usually used fo r SED fit- 
ting for high-z LAEs. E.g. lGuaita et al.l ((201 lh test three 
different star-formation histories when performing SED 
fitting to z ~ 2.1 stacked LAEs, and find equally good 
fits to the data. The authors also note that SED fitting 
can relatively well constrain stellar mass and dust extinc- 
tion among the free parameters. The derived properties 
are given in Table [4j and the best-fit SEDs are shown in 
Figure 

One problem is that COSMOS-30679 was found to be 
blended with another brighter adjacent object via exam- 
ination of the COSMOS HST/ACS F814W imaging data 
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Fig. 5. — Results of SED fitting. The filled squares show the observed flux densities used for the fitting (B, V, r, i, z, J, (H), K, [3.6], 
[4.5], [5.8], and [8.0]), while the open squares indicate those we omit for the fitting (U and NB387), due to the unknown flux contributions 
of IGM absorption. The red lines show the best-fit model spectra, and the red crosses correspond to the best-fit flux densities. 



TABLE 4 

Physical Properties from SED Fitting 



Object fss 


log(A^) 


E(B - V) 


log(age) 


log(SFR) 


x 2 




(M ) 


(mag) 


(yr) 


(M /yr) 





COSMOS-08501 


n n+ 10 
u - u -o.o 


7 SO+^" ^ 
' '° Z -0.2S 


0.08tg;S 






1.821 


COSMOS-13636 


U -'-0.5 


q qo+O.OS 


n 97 +0.02 
u - z '-0.08 


6.i8±i;^ 


! 9+1-16 
°' lz -1.28 


21.162 


COSMOS-30679 


1 0+°'° 


11 14+O.OO 


rs+o oo 

u ' D ' 3 -0.00 


6 42+ ' 00 


A 79+O.OO 

^•' z -o.oo 


2391.029 


COSMOS-30679(t) 


i o +0 ° 


10.301°;^ 


n ko+0.03 
u -°°-0.03 


6.42±°-°* 


oc+1.44 
3.»8_ 23 


19.769 


COSMOS-43982 


o 2+ 01 


ln 70+O.O8 


n 4Q+O.OO 
u - 4J -0.04 


ni+0.25 

°- J1 -0.00 


9 t-Q+0.00 


71.371 


CDFS-3865 


o 7+ 01 

u -'-0.1 


q Kn+0.03 
y - ou -0.02 


n 1S+ 001 
al8 -o.oi 


7 n9 +0.06 
'• uz -0.04 


9 4O+0.05 
z - 4y -0.06 


16.854 


CDFS-6482 


i n+ 00 


q 79 +0.09 
a -'^-0.07 


n 1s +0. 03 
U - i8 -0.02 


7 oe+0.20 
'• 8t5 -0.23 


1 qo+° 15 


4.986 


SSA22-8043 


o o +OA 
u - u -o.o 


in n7+°- 07 


03+ 002 

U - UJ -0.01 


q 44+0-00 
y -*^-0.10 


n 76+ - 07 

u -'°-0.03 


112.152 



Note. — Physical properties and their lo" errors of LAEs from SED fitting. Stellar mctallicity is fixed to 0.2 Zq . 
The second column shows the escape fraction of ionizing photons (a free parameter), (f) Values obtained 

after removing the contribution from an adjacent object ( §4^2) . 



COSMOS-08501 



COSMOS-13636 



COSMOS-30679 



C0SM0S-439B3 




Fig. 6.— HST/ACSF814W images for COSMOS objects. White 
tick indicates 1 arcsec, corresponding to ~ 8.3 kpc at 2 = 2.2. 
North is up and east is to the left. 

(Figure IB])- In order to deblend this object, we use the 
galaxy profile fitting software GALFIT (V3.0; |Peng et al.l 
120101) . We first run GALFIT in the B band and fit 
both objects simultaneously. Then, we perform fittings 
in other bands by fixing parameters (except for bright- 
nesses) to those derived in the B band fitting and derive 
the deblended photometry, which is listed in Table 
We confirm that the i' band brightness estimated from 
the deblended object is consistent with the photometry 
in the HST/ACS F814W high-resolution image within 
their errors. We add residuals after subtracting the fit- 
ted profiles into the photometric errors in quadrature, so 
that those errors are used in the SED fitting. Since the 
adjacent object appears dominant at longer wavelengths 
and deblending gets less reliable, we do not perform the 



deblending for the IRAC images. The best fit SEDs of 
COSMOS-30679 after and before deblending are shown 
in upper right two panels in Figure [5] Although the 
bumpy shape of the SED after the deblending suggests 
additional contamination from the adjacent object, the 
fitting is acceptable (relatively small \ 2 value) compared 
with that obtained before deblending. In the following 
sections, we use the quantities after deblending. 

We obtain from our sample stellar masses ranging from 
7 x 10 7 M Q to 2 x 10 10 M Q (except for the AGN LAE) and 
dust extinction E{B - V) = 0.03 to 0.53. These results 
are consistent with the previously recognized trend that 
LAEs at z ~ 2 are diverse, unlike higher redshift coun- 
terparts (e.g., iNilsson et al.l |20TTI ) . We also find that 
dust extinctions estimated from the Ha/H/3 ratio and 
SED fitting are consiste nt with each other for CDFS-3865 
(Hashimot o et al.l l2012). This supports our assumption 
that the dust extinction from SED fitting can be rea- 
sonably applied for the extinction correction of nebular 
lines. 



4.3. Ionization Parameter and Metallicity 
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TABLE 5 

Physical Properties from the Nebular lines 



Object 


12 + 
N2 


log(0/H) 
R23 


(10 7 emu- 1 ) 


SFR 

(M yr- 1 ) 


SFR C or 

(Mq yr" 1 ) 


EW(Ha) 

(A) 


L(Ha) 

(10 42 ergs- 1 ) 


EW(Lya) 

(A) 


L(Lycv) 

(10 42 cr gB -l) 


AGN 

(10) 




(i) 


(2) 


(3 


) 


(4) 


(5) 


(6) 


(7) 


(S) 


(9) 




COSMOS-08501 


< 8.73 










6 7+ 2 ' 7 

"■'-2.1 


> 280 


0.67 + 0.13 


255 ± 26 


8.41 + 0.40 





COSMOS- 13636 

COSMOS-30679(t) 

COSMOS-43982 


< 8.61 

< 8.70 


8.08-8.67 


3.4- 


15 


7.5li;j| 

Q 0+ 08 
a ' u -0.8 

19.01 J;? 


17 q+2.6 
A '-3-4.9 
45 2 +5 ' 4 
71 9+6.2 


go + 26 

no+25 
ao -33 

4ll 4 4 


0.95 + 0.13 
1.14 + 0.10 
2.41 + 0.21 


73 + 5 
87 + 7 
105 + 8 


11.35 + 0.50 
8.47 + 0.65 
11.00 + 0.47 





1 


CDFS-3865 
CDFS-6482 
SSA22-8043 


< 8.35 

< 8.72 

< 8.90 


7 84+ ' 24 


24j 


20 

-9 


W7.9± 6 6 - 7 7 

Qy tj+5.0 
z ' • -5.0 
41 0+ 120 


190.0li29 

A8.0+H 
44 5+13.2 


oiq+175 
2fil +103 

1 20+ 53 
lzu _69 


13.65 + 0.85 
3.45 + 0.63 
5.19 + 1.52 


64 ± 29 
75 ± 52 
28 + 4 


29.79 + 4.93 
15.40 ± 8.09 
2.22 + 0.38 








Note. — (1) Oxygen abundance derived from _/V2-indcx ( fl4.3t . If [Nil] is undetected, la upper-limit is given. (2) Oxygen abundance derived 
from i?23-index ( j]4.3t . Ionization parameter is derived at the same time. (3) Ionization parameter (10 7 cms _1 ) estimated from the [O III] /[Oil] ratio 
and Z(R23) ( A4.3H . (4) Star formation rate (Mq yr -1 ) derived from uncorrected Ha luminosity f jj4.4| >. (5) Star formation rate {Mq yt~ ) derived 
from Ha luminosity after correcting for extinction derived from SED fitting f CTIt . (6) Rest-frame equivalent width of Ha (A; i|4~5t . (7) Observed 
luminosity of Ha (10 42 ergs -1 ). (8) Rest-frame equivalent width of Lya (A; see Table (9) Observed luminosity of Lya (10 42 ergs -1 ). (10) Flag 

of AGN ( Ej4.H l: AGN candidates arc flagged with "1". (j") Values obtained after removing the contribution from an adjacent object f ^4.2t . Since H/3 
and [Oin]A4959 arc inferred from other lines and thus involve large uncertainties, its mctallicity (both iV2 and i?23) and ionization parameter are 
given with la uncertainties. 



The ionization state in an Hn region is often char- 
acterized by the ionization parameter, q^J^l which 
is the ratio of the mean ionizing photon flux to the 
mean hydrogen atom density. Larger q ion means the 
Hll region is more highly ionized, followed by strong 
emission lines with higher ionization potentials. In 
this sense, the ionization parameter is well determined 
using ratios of emission lines of different ionization 
stages of th e same element, s uch a s the [O in] /[On] 
ratio (e.g., IKewlev fc Dopital 12002ft . However, the 
[O in] /[On] ratio depends not only o n ionization state 
but a lso on gas metallicity (Fig. 1 of IKewlev fc Dopital 
l2002fh Conversely , metallicity estimates using i?23-index 
( Page! et alJll979l) : 

[On]A3727+ [Om]AA5007,4959 
R23= H£ (1) 

are a ffected by the ionization parameter to som e extent 
(e.g.. IKewlev fc DopitalfeOOa iNagao et al.ll2006D . There- 
fore, we combine the [Om]/[Oii] ratio and i?23-index 
(both corrected for dust extinction estimated from SED 
fitting) to estimate its ionization parameter and metal- 
licity at the same time iteratively. 

We describe the procedure for CDFS-3865 as an ex- 
ample. As an initial guess of metallicity, we use 
the empirical i?23 indicat or determined by local galax- 
ies (|Maiolino et al.l 120081) . We obtain two solutions 
12 + log(0/H) = 7.45±{J:ig and 8.66±°J°. We then 
use another altern ative indicator, jV^ -index defined 
as [Nn]A6584/Ha (jMaiolino et al.l 120081). The upper- 
limit of [Nil] /Ha provides an upper- limit of metallic- 
ity 12 + log(0/H) = 8.35 (la-), which removes the 
high-metallicity solution by i?23-index at the 2.5a level. 
Therefore, the empirical relation yields the metallicity 
12 + log(0/H) = 7.45±g:?|. As the second step, we use 
this metallicity and t he [Qui] /[On] ratio to es timate its 
ionization parameter. IKewlev fe Dopital (|2002| ) calculate 
the relations between [O in] /[On] ratio and ionization 

15 We use the subscript "ion" for the ionization parameter to 
distinguish it from the g-parameter which st ands fo r the effect of 
dust on Lyo compared with UV-continuum ( i]5,2.1[l . 



parameter with discrete values of metallicity. Among 
them, our initial guess of metall icity corresponds to 
their lowest metallicity (~ 0.09 Zq; lAllende Prieto et al.1 
2001). By using the relation, we find that the high 
[O in] /[On] ratio can only be reproduced with ioniza- 
tion parameters as high as several times 10 8 cms . 
As the third step, w e reca lculate its metallicity from 
the IKewlev fc Dopital (|2002t )'s relation (i?23-index) for 
the ionization parameter 3 x 10 8 cms~ 1 , and obtain 
12 + log(0/H) = 7.84l~°Jf3 Finally, we recalculate its 
ionizati on parameter based on t he metallicity and Eq. 
(12) of IKewlev fc Dopital (j2002l f 7 l. and obtain q ion = 
2.4j^Q g x 10 8 cms -1 . This is consistent with the assumed 
ionization parameter in the third step (3 x 10 8 cms _1 ). 
We check that further iterations do not change our final 
solutions. To summarize, CDFS-3865 has the ionization 
parameter qi on — 2.4l" 'g x 10 8 cms -1 and the metallicity 
12 + log(0/H) =7.84^. 

In the same manner, we estimate COSMOS-30679's 
ionization parameter and metallicity to be q.i on = (3.4— 
15) x 10 7 cms- 1 and 12 + log(0/H) = 8.08-8.67, re- 
spectively. Since COSMOS-30679 does not have H/3, its 
strength is inferred from Ha assuming the Case B and the 
dust extinction. In addition, its broadband SED is fitted 
after being deblended with the adjacent object ( §4.2| . 
Therefore, the object may possess larger uncertainties. 
To be conservative, we adopt the 2a range as the un- 
certainty for its derived quantities. We check that the 
upper-limit of metallicity from the A2-index is consis- 
tent within the 2a level. 

For the LAEs without [O II] , we estimate their metal- 
licities by using the empirical A2-index relation, and no 
constraint on ionization parameters is provided for them. 
The metallicities and the ionization parameters are sum- 
marized in Table [5] 

4.4. SFR 

10 Using the 1.5 X 10 8 cm s _1 relation does not change the metal- 
licity estimate significantly (~ 0.02 dex lower). 

17 Although the original relations are suspended at qi on = 3 X 
10 s cms -1 , we extend their polynomials toward higher qi on values. 
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The Ha luminosity is thought to be the most reliable 
SFR indicator relative to those based on the rest-frame 
UV and optical spectral features. Indeed, it is propor- 
tional to the birth rate of massive stars as well as being 
relatively insensitive to dust extinction as compared with 
UV-continuum. We measure the SFR of the LAEs fro m 
the Ha luminosity using the relation (|Kennicu 

SFRfMoyr" 1 ] = 7.9 x 10~ 42 L(Ha) erg s -1 . (2) 
In Table O we summarize the derived SFRs. 

4.5. Ha Equivalent Width 

We calculate the Ha equivalent width of an object from 
the Ha flux divided by its continuum flux density derived 
from K band photometry. To estimate the proper con- 
tinuum, however, we have to remove the contribution of 
the Ha emission on the K band photomet ry. We follow 
the pr escriptio n given in the appendix of iGuaita et al.l 



tne p rescription give n m tne a ppendix 
(|20TO) (see also IFinkelstein et alllMlj ): 



/„(Ha) =R T x 



F(Ha) 



/ Tk (A) j2 d\ ' 



(3) 



where /^(Ha) is the amount of flux that Ha line con- 
tributes to the K band photometry, F(Ha) is the ob- 
served Ha flux, Rt is the ratio between the filter trans- 
mission at Arq and the maximum of the K band trans- 
mission function, and Tk (A) is the K band filter trans- 
mission at a given wavelength. The flux estimated from 
the calculation is subtracted from the K band photom- 
etry The corrections span from 0.05 (COSMOS-43982) 
to as much as 0.63 mag (CDFS-3865). Since COSMOS- 
08501 is not detected in the K band, we use the la lower- 
limit of the K band photometry. We thus obtain an 
upper-limit of the continuum (correction > 0.26 mag), 
hence a lower-limit of the EW(Ha) for COSMOS-08501. 
The derived EWs(Ha) are listed in Table [5j 

5. DISCUSSION 

In this section, we discuss the physical properties of 
LAEs determine d by this study in conjunction with 
IFinkelstein et al.l ([20111 ) and lNakaiima et all (|2012t ). In 
the first subsection, we make comparisons in terms of ion- 
ization parameter, metallicity, and SFR of LAEs at z ~ 2 
as well as lower-redshifts, and other galaxies at similar 
redshifts such as LBGs. We then examine the implica- 
tions. In the second subsection, we show the properties 
of LAEs newly found by spectroscopy based on Lya and 
Ha hydrogen recombination lines. We extend the dis- 
cussion to understanding the origins of the strong Lya 
emission observed in LAEs. 

5.1. Comparisons between LAEs and other galaxies 
5.1.1. Ionization State 

We find high ionization parameters for the LAEs. In 
this section, we examine this point further by comparison 
to other galaxies, and discuss its implications. 

The situation is most clearly seen in Figure which 
shows the [O in] /[On] ratio versus i?23-index. We plot 
CDFS-3865 and COSMOS-30679 (red circles), and other 
galaxies such as LBGs at z = 2-3 (blue symbols). The 
black grid represents model predictions of [O in] /[Oil] 
ratio and i?23-index at a given metallicity and (discrete) 



ionization parameter ([Kewlev fc Dop ita 2002). Accord- 
ing to the calculation, CDFS-3865 is close to the qi on = 
3 x 10 s cm s _1 curve. This confirms its high ionization pa- 
rameter and low oxygen abundance. The high ionization 
parameter is suggestive of CDFS-3865 showing a very 
hard ionizing spectrum. CDFS-3865 is therefore consid- 
ered to be young galaxy dominated by massive stars in 
possibly small regions (expected for young Hn regions). 
Its low metallicity (Z ~ 0.1 Z e ) and large EW(Ha) 
(~ 800 A) also support this idea (see also 35.1.2p . 

Compared to CDFS-3865, most other z — 2-3 galax- 
ies appear to have lower ionizatio n parameters. Notabl e 
exceptions are o bjects found by IFosburv et al.1 ( 2003 ). 
| Erb et alj (I2010D. a nd ( partly) IRichard et al.l (|201lli . 
IFosburv et all (|2003l) and lErb et alj (|2010l) reveal from 
their multi-emission lines analysis that their galaxies pos- 
sess high ionization parameters (qi on ^ 10 9 cms~ 1 ), low 
metallicities (Z ~ 0.1 Z or less), being roughly com- 
parable to (or more or less extreme than) those found 
in CDFS-3865. Interestingly, both galaxies turn out to 
exhibit strong Lya emission. This supports our sugges- 
tion that strong L ya galaxies are in pa rt represented 
by young galaxies. IRichard et al.1 ([20111 ) also find two 
lensed galaxies (blue diamonds with upward arrow) pos- 
sibly have high ionization parameters (qi on ^ 10 s cms -1 ) 
judged from their metallicities determined empirically 
(12 + log(0/H) = 8.00t io and 8 - 77 -o.it for Sextet and 
MACS0712, respective ly). Unfortuna t ely, L ya measure- 
ments are not given in IRichard et all poll . 

However, as indicated by COSMOS-30679, LAEs' ion- 
ization parameters are unlikely to be always very high. 
The ionization parameter for COSMOS-30679 is esti- 
mated to be qi on = (3.4-15) x 10 7 cm s" 1 , comparable to 
those found in other high-z galaxies (e.g.. | Hainline et aT] 
[20091 ) or LBGs (e.g.. iMannucci et al.ll2009D . Therefore? 
the ionization parameter is not likely the quantity that 
distinguishes LAEs from non-LAEs. 

Next, we compare our LAEs with galaxies s een i n 
the local universe. According to iDopita et al.1 (2006), 
the SDSS galaxies have ionization parameters qi on = 
several x 10 7 cms~ 1 on average. High-z galaxies includ- 
ing the LAEs appear to have systematically higher ion- 
ization parameters than the S DSS galaxies. Ga laxies at 
intermediate redshift (z ~ 0.7; ILillv e t al. 2003|) are also 
plotted in Figure [7J They show systematically higher 
[Oiii]/[Oii] ratios than the SDSS galaxies. Since their 
metallicities are comparable to those of the SDSS galax- 
ies, the offset can be interpreted that higher- z galaxies 
tend to have higher ionization parameters. 

Another interesting comparison is with low-metallicity 
galaxies in the local universe (gray diamonds). Un- 
like typical SDSS galaxies, they appear in almost the 
same parameter space occupied by high-z galaxies on 
Figure \7\ Since their low-metallicities (7 < 12 + 
log(Q /H) < 8.5) are de termined by direct T e method 
(e.g., Izo tov et al.1 I2006D and thus reliable, the high 
[O in] /[On] ratios can be int erpreted as due to high 
ionization parameters. Indeed, iNagao et al.1 ([2006D sug- 
gest that low-metallicity galaxies have higher ioniza- 
tion parameters (> 3 x 10 8 cms _1 ) than more metal- 
rich galaxies seen in the SDSS sample. Therefore, 
these low-metallicity galaxies in the local universe ap- 
pear to have similar metallicities and ionization pa- 
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FlG. 7.— The emission line flux ratio [O m]/[0 il] vs. #23-index for our z = 2.2 LAEs (CDF S-3865: the solid red circle, COSMOS-30679: 
the open red circl e) and non-LAEs (th e blue symb ols); the metal-poor BX galaxy at z = 2 .3 HErb et al.ll2uTol : the open pentagon), lensed 
galaxies at z ~ 2 (Ha inlin e et aI.II 2009l: the crosses. iRichard et a.1 .11 20 111 : the open dia monds, Rigby ct al. 2 0111 : t he open square), the le nsed 
galaxy at z = 3.4 (Fosburv ct al. 2003: the solid star), z ~ 3 LBGs (Pcttini ct al. 2001: the solid pentagons, Mannucci ct al. 2009: the 
solid squares), as well as the star forming galaxies at 0.47 < z < 0.92 JLillv et al. 2003: the black open triangles), and local star forming 
galaxies in the low-metalli city range (7 < 12 + log(0/H) < 8.5: the gray diamonds) and high-metallicity range (12 + log( Q/H) > 8.2: 
the gray dots) compiled by Nagao ct al. (2006). The values on this plot are not corr ected for dust extinction , except for the lNagao et all 
(2006)'s local data points. The black grid shows photo-ionization model calculations (Kcwlcy & Dopita 2002). For reference, the relations 
of q ion = 3 X 10 8 , 8 X 10 7 , and 2 X 10 7 cms 1 are shown, linked to each other at same metallicities, which are denoted by 12 + log(0/H) 
values with dots. 



rameters as high-z star-forming galaxies, and particu- 
larly those with extreme quantities are likely analogs 
of high-z LAEs in terms of their high ionization pa- 
rameters and low metallicities. Among them, star- 
forming luminous compact galaxies (L CGs), or popularly 
referr ed to as "green pea" galaxies (jCardamone et afl 
2009) may indeed be good analogs. Although local 
low-metallicity galaxies tend to have much mo re modest 
SFRs than those observed for the LAEs (e.g., iLee et ail 
I2004D. LCGs have comparable SFRs similar to LAEs 
( Izotov et al.l 1201 If). Based on their low- metallicities 
(jAmorin et al.l 120101 ) . stron g [Oili] lines , and distribu- 
tions in the BPT-diagra m (jlzotov et al.l 120111 : see also 
Fig. 5 of lErb et al.l [2010 ) . LCGs are likely to have high 
ionization parameters (> 10 s cms -1 ). In the SDSS 
spectroscopic galaxy sample, 139 galaxies are in the 
low-metallicity ra nge while 48, 497 g alaxies in the high- 
metallicity range (jNaeao et al.l [2006^. Therefore, the 
local low-metallicity galaxies occupies only ~ 0.3 % of the 
SDSS galaxy sample. Similarly, LCG s occupies merely 
~ 0.0 7% of the SDSS galaxy sample (jCardamone et all 
2009). Based on t he luminosity functi on determined for 
the SDSS galaxies (Bl anton et al.ll2003h . the number den- 

18 Each sample has slightly different conditions such as S/N 
ratios of lines or redshifts. The direct comparison with the two 
samples' numbers is not entirely self-consistent, but provides an 
approximate ratio. 



sity is ~ 1.5 x 10 -2 Mpc~ 3 above a detection limit M r = 
— 17.63, which corresponds to 0.037 L*. The limit is cal- 
culated from the app arent magnitude limit (r < 17.77; 
iTremonti et al.l [20041) for the spectr oscopic sample an d 
the redshift threshold (z > 0.028; iNagao et all 120061) . 
The number density of the local low-metallicity galax- 
ies and LCGs are thus approximately 4.4 x 10~ 5 and 
1.0 x 10~ 5 Mpc~ 3 , respectively. On the other hand, the 
number density of LAEs at z = 2.2 is ~ 1.7 x 10~ 3 Mpc -3 
calculated by integra t ing t he luminosity function pro- 
vided bv lHaves et ail (UlloQ above L = 0.037 L*. The 
number density of the local low-metallicity galaxies is 
thus more than one order of magnitude smaller than that 
of LAEs. LCGs seem to be much less abundant by a fac- 
tor of about four. Such low-metallicity galaxies may be 
much more abundant at high-z: while rarely seen in the 
local universe. 

The final noteworthy implication of the high ioniza- 
tion parameters found in high-z LAEs is that such galax- 
ies may provide additional photons that cause hydrogen 
reionization of the IGM in the early universe. Previous 
censuses of early galaxies have revealed a possible short- 
age of ionizing ph o tons for the cosmic rei on ization (e.g., 
Ouchi "et~aT1 [200l iRobertson et al.l [2010h . lOuchi et al.1 
(2009) find that the universe could not be totally ionized 
by only galaxies at z = 7 if there is no evolution of prop- 

19 The detection limit is L = 2.8 X 10 41 ergs -1 (~ 0.02 L*). 
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Fig. 8. — (left) Mass-metallicity (M-Z) relation for z ~ 2 LAEs (red symbols), lensed galaxies, and continuum selected galaxies (blue). 
The circles show our LAEs, whose metallicities are measured/constrained by using i?23-index (with error bar)/7V2-index (lcr upper-limits). 
The open circle shows COSMOS-30679, whose metallicity range is given with 2<r. The red dot-dashed line s hows the la upper-li mit of 
metallicity (A r 2-index) for th e composite spectrum in K band. We additionally plot the average LAE (pentagon: Nakajima ct al. 2012) and 
two z ~ 2.4 LAEs (triangles: Finkclstcin ct al. 2011). LAEs at z = 0.195 — 0.44 are also show n with gray enclosed area {Cowic ct al .1 120111) . 
For non-L AEs, we plot len sed galaxies at 1.5 < z < 2.5 (open di amonds: [R ichard ct al. 20TJ), the metal-poor BX galax y at z = 2.3 (op en 
pentagon: Erb ct al. 2010), and BX/BM galaxies (open squares: |Erb ct al. 2006a|). The M-Z relation compiled bv lMaiolino et al. (2008) at 
z ~ and z ~ 2 are drawn as long-das hed and da shed curves, re spectively. All data plotted here have been recalibrated to have the same 
met allicity scale ( Maiolino ct al. 2008) and IMF {Salpctcr 1955) except for our LAEs with R 23-index whose metalli cities are calibrated 
with Kewley & Dopita (2002)'s relations, (right) Fundamental metallicity relation proposed by Mannucci ct al. (2 0TfT). The dashed curve 
and t he shaded areas indicate the relation and its typical dispersions respectivel y, determined in the lo cal u niverse ( Mannucc i et al.1120101 . 
120111) . The open squares show z ~ 2 star forming galaxies collected/compiled by Mannucci et al. ( 20101) from lErb et al.l H2006al ) and so on. 
Other symbols are the same as those in the M-Z relation. 

erties (e.g., escape fraction of ionizing photons, metallic- 
ity, dust extinction) from z ~ 3 to z ~ 7. Since the ion- 
ization parameter is by definition a measure of the num- 
ber of ionizing photons per hydrogen atom, the Hn re- 
gions are expected to be density-bounded when it is very 
high. In the conditions of density bounding, more ioniz- 
ing photons can escape from the Hii regions when the 
ioniz a tion parameter is highe r (e.g, iBrinchmann et aLI 
20081: iGiammanco et al.l [2005). Therefore, the escape 
fraction of ionizing photons is expected to be higher for 
galaxies (assembly of the Hn regions) with higher ioniza- 
tion parameters. Since the fraction of LAEs among star- 
forming galaxies i s known to increas e with redshift (e.g., 
lOuchi et all 120081 : iStark et al.l 120101 ) . LAEs and similar 
low-mass galaxies could produce ionizing photons effi- 
ciently and play a key role in supplying ionizing photons 
for the cosmic reionization at z > 6. 

5.1.2. Metallicity 

Figure [8] (left) shows the observed mass-metallicity 
(M-Z) relation. The red symbols represent LAEs at 
z ~ 2, while the blue symbols show the other galaxies 
such as LBGs at the similar redshifts. Although some 
of our LAEs have just weak upper-limits of metallicity 
due to the lack of [N ii]-detection (g3J), for CDFS-3865 
and COSMOS-30679, we obtain the metallicity estimates 
from i?23-index with the oxygen lines. We plot these 
i?23 metallicity estimates with larger circles in Figure [5J 
Both of them fall below the the conventional M-Z rela- 



tio n of z ~ 2 LBGs. At lea st one LAE (HPS194) found 
by iFinkelstein et all (|2011| ) also looks less chemically- 
enriched for its mass. Although the offsets of these 
LAEs from the M-Z relation can be due to the intrin- 
sic scatter of the rela tion seen in the local universe (e.g., 
iTremonti et~all[200l . it may also indicate that they are 
less chemically evolved for their stellar masses. The idea 
is consistent with the suggestion of their high ionization 
parameters. Since galaxies with high ionization parame- 
ters (e.g., BX418) tend to fall below the M-Z relation as 
welE3, young galaxies with high ionization parameters 
may not follow the M-Z relation defined by more evolved 
galaxies. Alternatively, differences in star-formation ac- 
tivity may cause the scatter. We discuss this point later. 

We also constrain an average metallicity from the K 
band composite spectrum (Figure |4]) by using the em- 
pirical [Nn]-index. The la (2a) upper-limit of metallic- 
ity is 12 + log(0/H) < 8.42 (8.66), which corresponds 
to Z < 0.54 (0.93) Z e . On the other hand, we in- 
dependently obtain an average lower-limit of metallic- 
ity for LAEs to be 12 + log(0/H) > 7.93 (7.63), or 
Z > 0.17 (0.09) Z Q at the la (2a) level, based on the 
[Oii]/(Ha+[Nn]) ratio whose fluxes are obtained by 
stacking 1.18 an d 2.09 /urn narrowban d images for more 
than 100 LAEs (jNakaiima et al.ll2012l ). LAEs thus typi- 



20 Galaxies wit h high ionization parameter ( Fosbur v et al.120031 : 
I Richard ct al. 2011]) are omitted due to their high redshifts (z > 
2.5). Since the M-Z relation is known to evolve with redshift, a 
direct comparison with the z ~ 2 LAEs is not appropriate. 
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cally have a metallicity 12 + log(0/H) = 7.93-8.42 (7.63- 
8.66) at the la (2a) level. The range is robust in the 
sense that the upper-limit is constrained by bright, mas- 
sive LAEs while the lower-limit by faint, low-mass LAEs. 
Super-solar metallicities and zero metallicities of Pop III 
are thus in general unlikely for LAEs in this epoch. The 
metallicity range also suggests that LAEs at z ~ 2 are 
rather less chemicall y enriched than thos e at z ~ 0.3 
( 12 + log(0/H) ~ 8.4^1; ICowie etlH l20ll ■ 

INakaiima et all (|2012ft find that LAEs fall typically 
above the M-Z relation below the stellar mass ~ 10 9 Mq. 
In contrast, the current study finds less chemically en- 
riched LAEs for their masses. This apparent incon- 
sistency may be due to the sampling of a large va- 
riety of evolutionary phases within the LAE popula- 
tion. Such v ariation has been indeed reported by other 
studies (e.g.. iNilsson et al.ll201l1; lOteo et all 12012ft . Al- 
ternatively, differences in star-formation activity may 
cause the inconsistency. To check the possibility, we 
plot th e LAEs on the funda mental metallicity relation 
(FMR; IMannucci et al.H2010l ). the relation between stel- 
lar mass, metallicity, and SFR. In Figure [8] (right), most 
of the galaxies at z = 0-2 including the average LAE 
(|Nakaiima et al.ll2012ft appear to be consistent with the 
same FMR determined by the SDSS galaxies within their 
errors. However, some LAEs are significantly below the 
relation. We can speculate that the FMR is not universal 
and may fail to reproduce the properties of very young 
galaxies. Clearly however, much more data is needed to 
test the idea statistically. 

5.1.3. Star- Formation Activity 

In order to examine the star-formation activity of 
LAEs, we plot the LAEs on the sSFR versus stellar 
mas s plane in Figur e [9l We also plot BX/BM galax- 
ies (|Erb et al.ll2006bft and sBzK galaxies (Hayash i et al.l 
120091 ). We note that this plot should be interpreted with 
cares since the spectroscopic data introduce limits in the 
sensitivity to low SFRs. 

BX/BM and sBzK galaxies follow a simple scaling re- 
lation between sSFR and the stellar mass, whose tight 
relation is r eferred to as the star formation main se- 
quence (e.g.. lDaddi et al.ll2007f ). Note that BX/BM and 
sBzK galaxies appear to have slightly different sequences; 
BX/BMs show lower sSFRs at a given mass. Compared 
to them, LAEs appear to follow the BX/BMs' main se- 
quence almost over the full mass range. CDFS-3865 and 
COSMOS-30679, which exhibit high ionization parame- 
ters and low metallicities, are not outliers on this dia- 
gram. This trend indicates that star-formation activities 
are well determined by their stellar mass, irrespective of 
the presence of Lya emission. 

5.2. Physical Properties inferred from Lya and Ha 
emission 

5.2.1. Correlation between EW(Lya) and EW(Hct) 

Figure [TU] shows equivalent widths of Lya and Ha for 
LAEs (hereafter referred to as the "EWs diagram" ) . This 
plot is useful to understand the star-formation histories 
of LAEs, because the EWs' continuum fluxes evolve in 

21 We recalculate the metallicity by using the IMaiolino et al.l 
(2008) indicator. The original estimate is ~ 0.15 dex lower. 
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Fig. 9. — Relation between stellar mass and specific star for- 
mation rate (sSFR) for z ^ 2 LAEs (red) and continuum se- 
lected galaxies (b lue); LAEs with NIR sp ectroscopy (circles: this 
work, triangles: Finkclstcin et al. 20111). the a verage LAE ob- 
taine d by stacking NIR NB imaging (pe ntagon: I Nakaiima et al.l 
[2012]) , BX/BM galaxies (open squares: lErb et alj I2006bl) , and 
sBzK galaxies (solid squares: Hayashi ct al. 2009). COSMOS- 
43982, a possible object with AGN activity, is marked with black 
diamond. The red open circle denotes COSMOS-30679, which has 
less accurate SED fitting results due to the neighbor's contamina- 
tion. All sSFRs plotted here are derived from Ha luminosity with 
Salpeter IMF. The short horizontal bars in each symbols for the 
LAEs indicate the sSFR with no dust extinction correction. The 
dashed lines correspond to constant SFRs of 0.1, 1.0, 10, 100, and 
1000 Mgyr 1 . 

different ways when different star-formation histories are 
assumed. Another advantage of using the EWs is that 
both are pure observables and to zeroth order indepen- 
dent of dust extinction (the effect of dust extinction will 
be discussed at the end of this section). 

The superposed curves on Figure [TUl illustrate the evo- 
lutions of the EWs for the two extreme star-formation 
histories, instantaneous burst (dashed) and constant 
star-formation (solid) at several metallicities (jSchaererl 
20030- For the instantaneous burst, since very mas- 
sive stars (M* > 10 M@) complete their evolution within 
< lOMyr, both EWs decline rapidly. As a result, their 
curves evolve quickly to lower-left on the EWs diagram. 
For the constant star-formation, on the other hand, the 
EW(Lya) stops declining around ~ 100 Myr, because 
massive stars that are responsible for both Lya emis- 
sion and UV-continuum reach a steady mode. Since the 
EW(Hq) keeps declining as the older stars build up in 
the galaxy, their slopes of the evolutionary tracks be come 
less steep on the EWs diagram. Although iSchaeren s cal- 
culations stop at ~ 400 Myr for constant star-formation, 
the E W(Lyq) varies little after 100 Myr (|Charlot fc Falll 
11993ft . Therefore, the tracks must extend left almost hor- 

22 The models are collected from the Strasbourg astronomical 
Data Center (CDS). We present here t he three metallicity cases 
assuming Salpeter IMF (Salpeter 1955|) with upper (lower) mass 
cut-off to be 100 Mq (1 Mq). Case B recombination is assumed for 
an electron temperature of T e = 3 X 10 4 K at zero metallicity and 
T e = 10 4 K otherwise, and an electron density of n e = 10 2 cm" 3 . 
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Fig. 10. — Relation between Lya and Ha equivalent widths. The 
red symbols show z ~ 2 LAEs (s ame as in Figure [3, and the gray 
squares show LAEs at z = 0.3 (Cowic ct al. 2011). Superposed 
lines are evolutions of the EWs for instantaneous burst (dashed) 
and constant star-formation (solid) at metallicities of solar (or- 
ange) , sub-solar (green), and zero (blue) calculated by ISchaererl 
(2003). Ages are denoted by the numbers (6, 7, 8) near the dots 
on the lines which indicate (1 Myr, lOMyr, 100 Myr). 

izontally after the terminal points. 

Compared with the model predictions, some LAEs 

with EW(Ha) ~ 1000 A, CDFS-3865, HPS256 

(|Finkelstein et all l20ll . and possibly COSMOS- 
08501, can be explained by the instantaneous burst 
models with very young ages (a few Myr). The other 
LAEs whose EW(Ha) is modest are better described 
by the constant star-formation models with > 100 Myr. 
Although the sample is small, we find more than half 
of the LAEs are more consistent with a continuous 
star-formation history rather than an instantaneous 
burst. This finding supports our assumption of adopting 
a constant star-formation history in the SED fitting. 

Note, however, that since the EWs are sensitive to re- 
cent starbursts, LAEs on instantaneous tracks are not 
necessarily very young, but can be experiencing a burst 
after a continuous star-formation. Such a combination 
of burst plus continuous star-formation is indeed needed 
to explain the EWs of Ha and optical colors observed 
in local dwarf galaxies ([Led 120061 ). An LAE with such a 
combined star-formation history is on a track of instan- 
taneous burst during starburst phases and on a track 
of constant star-formation for the remaining time. The 
trend seen in the EWs diagram may be a result that most 
of the LAEs are not in the burst mode. 

Similarly, most of z ~ 0.3 LAEs appear to prefer con- 
tinuous star-formation on the EWs diagram. Interest- 
ingly, z ~ 2 LAEs appear to have systematically higher 
EW(Lya) and EW(Ha) than lower- z LAEs. Although 
our spectroscopic sample may be biased toward larger 
EW(Lya), the difference may be a sign that higher- z 
LAEs are younger. This idea is supported by the infer- 
ence of lower metallicities at higher- z ( H5.1.2j) . 

Although adopting an extended star-formation history 



complicates the estimates of metallicity and age from the 
EWs diagram, independent metallicity estimates would 
be helpful for constraining its star-formation history and 
age. E.g., we find CDFS-3865 has a sub-solar metallicity 
f £|4.3[) . The sub-solar, instantaneous burst track repro- 
duces its EWs well with an age of a few Myr. This very 
young age is consistent with what we infer for CDFS- 
3865 from its high ionization parameter and low metal- 
licity f £15.1. H and £J5.1.2[) . The age roughly agrees with the 
estimate from the SED fitting. Although some LAEs in- 
cluding COSMOS-30679 have inconsistent age estimates 
from the EWs diagram and the SED fitting, they can 
be due to the fact that ages are weakly c onstrained and 
lumin osity weighted by SED fitting ('e.g.. iGawiser et all 
l2006h . Particularly, COSMOS-30679 is thought to have 
much larger photometric errors ( £14.21) . Alternatively, 
their EW(Lya) may be affected by dust. We discuss 
the possibility further below. In any case, since age es- 
timates are only possible through SED fitting, the EWs 
diagram can provide another independent constraint on 
the stellar populations. 

One worry is the effect of dust on Lya. If the degree of 
dust extinction is different for Lya and UV-continuum, 
EW(Lya) is no longer independent of dust extinction. In 
order to examine this effect , we in troduce a q parameter 
following Finkclst ein et al.l ([2008D . The q parameter is 
defined as q — r(Lya)/ri2i6, where r(Lya) and T1216 are 
optical depth for Lya and UV-continuum at A = 1216 A, 
respectively. Small q values (< 1) mean Lya photons suf- 
fer less attenuation by dust than UV-continuum photons, 
while large values mean Lya photons are more heavily 
attenuated. In the former (latter) case data points go 
down (up) on the EWs diagram after the corrections. 
From pr evious works, LAEs at z ~ 2 have modest q val- 
ues; e.g.. lHaves et al.1 ( |2010P ) obtain q ~ 1-1.5 for z = 2.2 
LAEs whose Lya and Ha lum inosities are estimated from 
two narrowbands. Similarly, Nakaiima et al.l (|2012D ob- 
tain q = 0.7 ± 0.1. iBlanc et al l (|2011l) obtain q = 0.99 
for z = 2-4 LAEs whose intrinsic Lya luminosities are 
inferred from UV-continuum. Therefore, assuming q = 1 
is reasonable, though we should take care of the possible 
effect of dust on the EWs diagram. 

Another concern involves the validity of the Case B 
recombination assumption for the LAEs. We discuss the 
issue in the next section. Along with the discussion so 
far, however, the EWs diagram shows that LAEs can 
be well described by the conventional stellar population 
models and that there is no need to invoke non-standard 
IMFs. 

5.2.2. super Case B objects 

Figure [IT] (left) shows the Lya/Ha ratio against 
EW(Lya) for LAEs. The dark red symbols indicate the 
observed Lya luminosity divided by the intrinsic Ha lu- 
minosity corrected using the attenuation inferred from 
the SED fitting (i.e., showing lower-limits on the y-axis), 
while the light red symbols show the observed Lya/Ha 
ratios. A trend that LAEs with larger EW(Lya) have 
larger Lya/Ha ratio appears to be present. This trend 
itself is not so surprising, but interestingly some LAEs, 
COSMOS-08501, and possibly COSMOS-13636 and one 
LAE at z ~ 0.3, may have Lya/Ha ratios exceeding th e 
Case B recombination value (8.7; iBrocklehurstJ [19711) . 
which we call "super Case B". Super Case B Lya/Ha 
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Fig. 11. — (left) Lya/Ha flux ratio vs. EW(Lya). Symbols are the same as in Figure IT01 The darker symbols show the ratios calculated 
by the observed Lya luminosity divided by the intrinsic Ha luminosity (i.e., lower-limit of the ratios), while the lighter show the ratios 
calculated by the observed Lya and Ha luminosities. The ratios for z ~ 0.3 LAEs are based on observed luminosities. Dotted curve 
represents a single power-law fit given in Eq. (El . The dashed line and the gray shaded area show the Lya/Ha flux ratio assuming the Case 
B recombination (8.7; Brocklchurst 1971) and its variation with electron density (8.1-9.2 with n e = 10 2 -10 10 cm - 3 ; [Hummer fc Storey! 
119871) when T e = 10 4 K, respectively, ( right) Escape fraction of Lya photons (assuming the Case B) vs. dust extinct ion. The red squares 
show z = 2.2 LAEs (Hayes ct al. 2010), and the faint-colored red error bars show z = 2—4 LAEs (Blanc ct al. 2011). Superposed curves 
show the relations at a given q parameter (q = 0.0,0.1,0.5, 1.0,2.0, 10.0). 



ratios have also be en reported for some strong LAEs in 
the local universe (lHaves et al.1 [20071 lAtek et al.1 [200l 
lOti-Floranes et all 120121 and references therein) . The 
variation of the intrinsic Lya/Ha ratio (the gray shaded 
area) does not seem to be a significant issue. Although 
it is not obvious that super Case B objects are included 
in our sample due to the large errors, we discuss possible 
physical origins of them in case they really exist. 

We consider the possible effect of geo metry, and the 
kinematics of dust and gas in the ISM. INeufeldl (|l99l 
propose a clumpy, multi-phase ISM where gas and dust 
are gathered in clouds within a low-density medium. 
With such circumstances, Lya photons can be scat- 
tered at the surfaces of the clouds due to the resonant 
nature, while continuum photons would penetrate the 
clouds deeply. Since dust is contained in the clouds, 
Lya photons would have a much smaller chance of en- 
countering dust than other photons. In this scenario, 
a large EW(Lya) and Lya/Ha ratio can be observed. 
Alternatively, an outflow of the ISM can be a cause of 
strong Lya emiss ion Je.g-. lKunth et al.lll998l ). However, 
IHashimoto et al.1 ([2012) find an anti-correlation between 
EW(Lya) and Lya velocity offset for LAEs. Thus, out- 
flows do not appear to be a major mechanism for pro- 
ducing large EW(Lya) and Lya/Ha ratio. 

In order to clarify the effects of a potentially clumpy 
geometry of ISM on Lya/Ha ratios, we plot in Figure 
[TT1 (right) the relation between E(B — V) and the es- 
cape fraction of Lya photons (/ c L s f) under the Case B 
recombination assumption. We estimate fiZ a as 



show the relations at a given q parameter ( tj5.2.1l) ; 
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where subscripts 'int' and 'obs' refer to the intrinsic and 
observed quantities, respectively. The superposed lines 



where fci2is is an extinctio n coefficient at A = 1216A 
(11.98; ICalzetti et aT]|2000D . The clumpy geometry of 
ISM (or outflow) is favored by objects with q — 0-1. 
From Figure ITT1 (right), most of the LAEs presented here 
are located in the range q = 0-1; e.g., COSMOS-30679 
and COSMOS-13636 have large f^/ c a in spite of their 
moderate amounts of dust. SSA22-8043 has a large q 
parameter of ~ 10, and can be an exception if its Lya is 
heavily resonant-scatted by neutral hydrogen gas. 

A notable object is COSMOS-08501. It appears to fall 
above the q — line, where the clumpy ISM model does 
not work assuming Case B recombination, albeit with rel- 
atively large errors. Moreover, since it is inferred to pos- 
sess relatively small amount of dust, the large EW(Lya) 
owing to the clumpy ISM is unlikely. In case the ob- 
ject is really super Case B, the only remaining explana- 
tion is the Lya enhancement caused by collisional exci- 
tations. Due to the decreasing collisional strengths with 
increasing principle quantum number, collisional excita- 
tions can lead to Lya/Ha ratios over the Case B value 
(see also lOsterbrocld[l989f l. Shocks caused by interac- 
tions with other sources, AGN activity, supernova explo- 
sions, strong outflows or infall are possible candidates for 
the collisional excitations. 

Based on the HST images (FigureEJ), COSMOS-08501, 
a super Case B candidate, looks very compact and shows 
no sign of interactions. COSMOS-13636, which has a 
very small q parameter, shows two faint sources nearby 
within ~ 5kpc (projected) from the object. Its strong 
Lya emission can be (partl y) due to shocks c aused by in- 
teractions. As discussed bvlMori "et~aT1((200l . supernova 
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explosions can cause strong shocks, resulting in strong 
Lya emission. Although the authors intend to explain 
extended Lya blobs (~ 100 kpc) with high Lya luminosi- 
ties (~ 10 43 ergs _1 ), their basic ideas can be applied to 
normal LAEs . AGN activity and ou tflows seem unlikely 
(see 347fl and Hashimoto et al.ll2012t ). 

When shocks (even partly) contribute to emission lines 
of LAEs, estimates of physical quantities such as SFR 
(from Ha), metallicity (from A2-index), ionization pa- 
rameter (from [O ill] /[Oil] ratio), and dust extinction 
(from Balmer decrement) become less accurate. In par- 
ticular, the intrinsic Ha/H/3 value becomes larger when 
shocks are present, resulting in an overestimate of the 
abundance of dust. This effect may a lso help explain the 
prese nce of super Case B objects (e.g.. lOti-Floranes et al.l 
I2012D . Thus, the results presented in this paper may re- 
quire some corrections for the presence of shocks. We 
plan to address this in future work, through the simulta- 
neous application of photo-ionization and shock models, 
with deeper spectroscopy which will detect the weaker 
lines not fully detected in the observations presented here 
(e.g., Balmer lines, [On], [Om], [Nil], [S11]). 

A final remark is that from our discussions so far it 
is evident that Lya is not a robust indicator of SFR for 
LAEs. The observed data (z = 0-2) in Figure [TT] (left) 
are relatively well represented by a single power-law fit 

log(Lya/Ha) obs = (-1.08 ±0.29) 

+(0.87 ±0.19) x logEW(Lya),(6) 

which is shown by the dotted curve in Figure [TT] (left). 
According to this simple relation, SFRs of LAEs with 
EW(Lya) ~ 20 A (a typical threshold in narrowband 
searches) can be underestimated by a factor of about 10. 
Thus, SFRs estimated from Lya may involve a factor of 
~ 10 uncertainties intrinsically. 

6. SUMMARY 

We have presented NIRSPEC and MMIRS rest-frame 
optical spectra of seven Lya emitters (LAEs) at z = 
2.2, which are selected from our Subaru/Suprime-Cam 
NB387 survey in COSMOS, Chandra Deep Field South, 
and SSA22. Our first NIR spectroscopic result dis- 
cusses the kinematics of LAEs and is presented in 
lHashimoto et al.l (|l012). As a companion study, this pa- 
per presents mainly the ionization and chemical proper- 
ties of LAEs based on multiple nebular lines. Our sam- 
ple includes one possibly AGN-dominated galaxy, and 
six star-forming galaxies. Ha is detected in all six star- 
forming LAEs, while [Nn]A6584 is only detected in the 
galaxy with signs of AGN activity. Among the six star- 
forming galaxies, one (CDFS-3865) also has detections 
of [On]A3727, H/3, and [O m]AA5007, 4959, and another 
(COSMOS-30679) has detections of [0 11] and [O in] . Our 
deep J band spectroscopic observations provide the first 
[0 11] -detections for two individual LAEs at high-z. Our 
main results are summarized as follows. 

• The [Oiii]/[Oii] ratio vs. i?23-index diagram re- 
veals that CDFS-3865 has a very high ionization 
parameter (qi on = 2.4+q'q x 10 8 cms _1 ) and a 
low oxygen abundance (metallicity; 12 + log(0/H) 
= 7.84^025) in contrast with moderate values of 



other high-z galaxies such as LBGs. COSMOS- 
30679 also has a relatively high ionization parame- 
ter (qion = (3.4-15) x 10 7 cms _1 ) and a low metal- 
licity (12 + log(0/H) = 8.08-8.67). LAEs would 
therefore represent an early stage of galaxy forma- 
tion dominated by massive stars in compact star- 
forming regions. 

• Local low-metallicity galaxies show similar ioniza- 
tion parameters and metallicities to LAEs, and can 
be local analogs of high-z LAEs. Among them, 
star-forming luminous compact galaxies (LCGs) 
may indeed be good analogs, having comparable 
SFRs similar to LAEs. The local low-metallicity 
galaxies and LCGs occupies only ~ 0.3% and 
0.07% of the SDSS galaxy sample, respectively, 
and are much rarer by more than one order of mag- 
nitude than high-z LAEs. 

• Since the escape fraction of ionizing photons is ex- 
pected to be higher with higher ionization param- 
eters, LAEs and similar low-mass galaxies could 
produce ionizing photons efficiently and play a key 
role in supplying ionizing photons for cosmic reion- 
ization at z > 6. 

• Both CDFS-3865 and COSMOS-30679 fall below 
the mass-metallicity relation of LBGs at similar 
redshifts. Their low metallicities do not seem to be 
explained by their star formation rates being taken 
into account, albeit with their relatively large error. 
Interestingly, galaxies with high ionization param- 
eters tend to fall below the relation. Very young 
galaxies may not follow the relation determined by 
more evolved galaxies. 

• The composite spectrum independently provides 
an upper-limit on the metallicity of 12 + log(0/H) 
< 8.42 (< 8.66) at the la (2a) l evel. Combined 
with an lower-limit of metallicity (Nak aTkna et al.1 
l2012f >. LAEs typically have metallicities 12 + 
log(0/H) = 7.93-8.42' (7.63-8.66) at the la (2a) 
level. 

• In contrast to the large differences in ionization pa- 
rameters and metallicity between LAEs and LBGs, 
we find LAEs have similar specific star formation 
rates as BX/BM galaxies at a given stellar mass. 

• We show that the EW(Lya) vs. EW(Ha) dia- 
gram is useful for understanding the star-formation 
histories of LAEs, being relatively independent of 
SED fitting. According to the EWs diagram, more 
than half of the LAEs appear to need an extended 
star-formation such as a burst superimposed upon 
a continuous star-formation (exponentially declin- 
ing or constant) rather than the instantaneous 
burst alone. 

• LAEs with low q- values (q — 0-1) can be explained 
by the clumpy geometry of ISM. Interestingly, our 
sample may include objects with further enhanced 
Lya, which we call super Case B. If they really 
exist, the only possible explanation is the colli- 
sional excitations of Lya. Interactions with other 
sources and/or supernova explosions are possible 
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key events that may cause shock-induced collisional 
excitation. If such shocks play a role in enhance- 
ment of the Lya flux, physical quantities such as 
SFR, metallicity, ionization parameter, and dust 
extinction should be re-computed using a combina- 
tion of photo-ionization and shock-excitation mod- 
els. We plan to investigate the role of shocks fur- 
ther in future works. 
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